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We present searches for the rare decay modes — > e^e~ , — > ^.^ ii~ , and — >■ e*/i^ in 
continuum e"'"e~ — >■ cc events recorded by the BABAR detector in a data sample that corresponds to 
an integrated luminosity of 468 fb^^. These decays are highly GIM suppressed but may be enhanced 
in several extensions of the Standard Model. Our observed event yields are consistent with the 
expected backgrounds. An excess is seen in the — >■ fj-'^ channel, although the observed yield 
is consistent with an upward background fluctuation at the 5% level. Using the Feldman-Cousins 
method, we set the following 90% confidence level intervals on the branching fractions: B{D'^ — ^ 
e+e" ) < 1.7 X 10~^ , ^ ) within [0.6, 8.1] x lO"'^ , and ^ e^^^F ) < 3.3 x 10"^ . 



PACS numbers: 13.20.Fc,11.30.Hv,12.15.Mm,12.60.-i 
I. INTRODUCTION 



In the Standard Model (SM) , the flavor-changing neu- 
tral current (FCNC) decays — >• are strongly sup- 
pressed by the Glashow-Iliopoulos-Maiani (GIM) mech- 
anism. Long-distance processes bring the predicted 
branching fractions up to the order of 10^^"^ and 10"^'^ 
for e+e^ and D'^ M^/^^ decays, respectively 

These predictions arc well below current experimental 
sensitivities. The lepton-flavor violating (LFV) decay 
— 7> e^/i^ is forbidden in the SM. Several extensions 
of the SM predict — >■ branching fractions that 

are enhanced by several orders of magnitude compared 
with the SM expectations The connection between 
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t^l and D'^ — D'^ mixing in new physics models 
has also been emphasized ^] . 

We search for D'^ — )■ decays using approximately 
468 fb^^ of data produced by the PEP-II asymmetric- 
energy e^e~ collider 3] and recorded by the TiABAR. de- 
tector. The center-of-mass energy of the machine was at, 
or 40 MeV below, the T(4S') resonance for this dataset. 
The BaBAR detector is described in detail elsewhere 
We give a brief summary of the main features below. 

The trajectories and decay vertices of long-lived 
hadrons are reconstructed with a 5-layer, double-sided 
silicon strip detector (SVT) and a 40-layer drift chamber 
(DCH), which are inside a 1.5 T solenoidal magnetic field. 
Specific ionization {dE/dx) measurements are made by 
both the SVT and the DCH. The velocities of charged 
particles are inferred from the measured Cherenkov angle 
of radiation emitted within fused silica bars, located out- 
side the tracking volume and detected by an array of pho- 
totubes (DIRC). The dE/dx and Cherenkov angle mea- 
surements are used in particle identification. Photon and 
electron energy, and photon position, are measured by a 
CsI(Tl) crystal calorimeter (EMC). The steel of the flux 
return for the solenoidal magnet is instrumented with lay- 
ers of either resistive plate chambers or limited streamer 
tubes 0, which are used to identify muons (IFR). 
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II. EVENT RECONSTRUCTION AND 
SELECTION 

We form D" candidates by combining pairs of op- 
positely charged tracks and consider the foUowing final 
states: e+e^ , /i^/i^ , e^/i^ , tt+tt^ , and K^tt^ . We use 
the measured D° — >■ tt+tt" yield and the known — > 
7r~'"7r~ branching fraction to normalize our — 
branching fractions. We also use the tt+tt^ can- 

didates, as well as the — >■ i^~7r+ candidates, to mea- 
sure the probability of misidentifying a tt as either a /i 
or an e. Combinatorial background is reduced by re- 
quiring that the candidate originate from the decay 
D*{2QIQ)+ D° 7r+ We select D° candidates pro- 
duced in continuum e+e^ — )■ cc events by requiring that 
the momentum of the 13° candidate be above 2.4 GeV 
in the center-of-mass (CM) frame, which is close to the 
kinematic limit for B D*7r, D*+ D°Tr+. This re- 
duces the combinatorial background from e"'"e~ — >■ BB 
events. 

Backgrounds are estimated directly from data control 
samples. Signal Z?" candidates with a reconstructed 13° 
mass above 1.9 GeV consist of random combinations of 
tracks. We use a sideband region above the signal re- 
gion in the mass ([1.90, 2.05] GeV) in a wide Am = 
m{D° n+) - m{D") window ([0.141, 0.149] GeV) to esti- 
mate the amount of combinatorial background. The 
and Am mass resolutions, measured in the D'^ — tt+tt^ 
sample, are 8.1 MeV and 0.2 MeV, respectively. We esti- 
mate the number of D'^ — >■ 7r"'"7r" background events se- 
lected as D'^ — > £^£^ candidates by scaling the observed 
D'^ — >■ 7r"'"7r~ yield, with no particle identification criteria 
applied, by the product of pion misidentification proba- 
bilities and a misidentification correlation factor G. The 
misidentification correlation factor G is estimated with 
the D'^ — >■ _ftr~7r+ data control sample. 

The tracks for the 1)° candidates must have momenta 
greater than 0.1 GeV and have at least 6 hits in the 
SVT. The slow pion track from the D*+ D" n+ de- 
cay must have at least 12 position measurements in the 
DCH. A fit of the D*+ D° n+-,D° ^ t+t' decay chain 
is performed where the tracks [t) are constrained 
to come from a common vertex and the D° and slow 
pion are constrained to form a common vertex within 
the beam interaction region. The probabilities of 
the D" and D* vertices from this fit must be at least 
1%. The reconstructed mass m{D^) must be within 
[1.65,2.05] GeV and the mass difference Am must be 
within [0.141,0.149] GeV. We subtract a data-Monte- 
Carlo difference of 0.91 ± 0.06 MeV, measured in the 
TT+TT^ sample, from the reconstructed mass 
in the simulation. 

We use an error-correcting output code (ECOC) al- 
gorithm jyj with 36 input variables to identify electrons 
and pions. The ECOC combines multiple bootstrap ag- 
gregated Q] decision tree @ binary classifiers trained to 
separate e, tt, and p. The most important inputs for 
electron identification are the EMC energy divided by 



the track momentum, several EMC shower shape vari- 
ables, and the deviation from the expected value divided 
by the measurement uncertainty for the Cherenkov an- 
gle and dE/dx for the e,7r,K, and p hypotheses. For 
tracks with momentum greater than 0.5 GeV, the elec- 
tron identification has an efficiency of 95% for electrons 
and a pion misidentification probability of less than 0.2%. 
Neutral clusters in the EMC that are consistent with 
Bremsstrahlung radiation are used to correct the momen- 
tum and energy of electron candidates. The efficiency of 
the pion identification is above 90% for pions, with a kaon 
misidentification probability below 10%. 

Muons are identified using a bootstrap aggregated de- 
cision tree algorithm with 30 input variables. Of these, 
the most important are the number and positions of the 
hits in the IFR, the difference between the measured and 
expected DCH dE/dx for the muon hypothesis, and the 
energy deposited in the EMC. For tracks with momen- 
tum greater than 1 GeV, the muon identification has an 
efficiency of around 60% for muons, with a pion misiden- 
tification probability of between 0.5% and 1.5%. 

The reconstruction efficiencies for the different chan- 
nels after the above particle identification requirements 
are about 18% for e+e", 9% for 13% for e^fi^, 

and 26% for tt+tt". The background candidates that 
remain are either random combinations of two leptons 
(combinatorial background), or 13° — > tt+tt" decays 
where both pions pass the lepton identification criteria 
(peaking background). The — > tt+tt" background is 
most important for the 13° — ;> /^^M^ channel. 

Figure [T] shows the reconstructed invariant mass dis- 
tributions from Monte Carlo (MC) simulated samples 
for the three D° — >■ £~^£~ signal channels. Also shown 
are the distributions from 13° tt+tt" reconstructed 
as — > and D° — ^ K^tt^ reconstructed as 

Z3° — > £^£^ for each signal channel. The overlap between 
the 13° — > and 13° — >■ tt+tt^ distributions is largest 

for the Z3° — > fX^ ii~ channel, while the 13° and 
Z3° K^TT^ distributions are well separated. 

The combinatorial background originates mostly from 
events with two semileptonic B and/or D decays. The 
sample of events selected by the above criteria are domi- 
nantly from e'^e^ — > BB events, rather than events from 
the e"^e~ — >■ qq, {q ~ u, d, s, c) continuum. We use a lin- 
ear combination (Fisher discriminant [lo[ ) of the follow- 
ing five variables to reduce the combinatorial BB back- 
ground: 

• The measured 13° flight length divided by its un- 
certainty. 

• The value of | cos 6'hoi | , where 6'hoi is defined as the 
angle between the momentum of the positively- 
charged 13° daughter and the boost direction from 
the lab frame to the Z3° rest frame, all in the Z3° 
rest frame. 

• The missing transverse momentum with respect to 
the beam axis. 
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FIG. 1: Reconstructed D^' mass (left) and Am (right) for the 
three signal channels: e^e~ (top), — >■ fj,'^ fJ-~ (mid- 

dle), and — >■ e^fj,^ (bottom). The solid (black) histogram 
is the signal MC, the dashed (blue) histogram is — >■ tt^tt" 
MC reconstructed as D° t~ , and the dotted (red) his- 

togram is D° — >■ K'-k'^ MC reconstructed as D° ^> l~ . 
The ^ tr and D° ^ TT+TT distributions have been 
normalized to unit area. The — > K~-k^ normalization is 
arbitrary. 



• The ratio of the 2'^'^ and 0*'' Fox- Wolfram mo- 
ments [iTI |. 

• The momentum in the CM frame. 

The flight length for combinatorial background is sym- 
metric about zero, while the signal has an exponential 
distribution. The |cos0hei| distribution is uniform for 
signal but peaks at zero for combinatorial BB back- 
ground. The neutrinos from the semileptonic decays in 
BB background events create missing transverse momen- 
tum, while there is none for signal events. The ratio of 
Fox- Wolfram moments uses general event-shape informa- 
tion to separate BB and continuum qq events. Finally, 
the signal has a broad 13° CM momentum spectrum that 
peaks at around 3 GeV, while combinatorial background 
peaks at the minimum allowed value of 2.4 GeV. 

Figure [5] shows distributions of the Fisher discriminant 
(T) for samples of BB MC, D° ^ signal MC, 

and continuum background MC. The separation between 
signal and BB background distributions is large, while 
the signal and continuum background distributions are 




Fisher Discriminant 



FIG. 2: Fisher discriminant, T , distributions for samples of 
BB MC (dashed blue), signal MC (solid black), 

and continuum MC (dotted red). The T distributions for 



— e^e and — > e*y^ are similar to those of 



TABLE I: Selection criteria for the three D° -s> l~ signal 
decay modes. The parameter in the last row is defined as 
(5Am = Am — Amo, where Amo is the nominal _D*+ — 
mass difference 



Parameter 






e /i+ 


cos 6'hci 


< 0.85 


< 0.90 


< 0.85 


T 


> 0.00 


> -0.25 


> 0.00 


m,(BP) (GeV) 


[1.815, 1.890] 


[1.855, 1.890] 


[1.845, 1.890] 


|(5Am| (MeV) 


< 0.5 


< 0.5 


< 0.4 



similar. For example, requiring T to be greater than 
removes about 90% of the BB background while keeping 
85% of the signal. The minimum J- value is optimized 
for each signal channel as described below. 

We use the | cos 6'hci | variable directly to remove con- 
tinuum combinatorial background. Figure |3] shows dis- 
tributions of I cos ^hci I before making a minimum !F re- 
quirement, for BB background, continuum background, 
and signal. The drop-off for | cos^hcil near 1.0 in the sig- 
nal distributions is caused by the selection and particle 
identification requirements. The BB background peaks 
near zero, while the continuum background peaks sharply 
near one. 

The selection criteria for each signal channel were cho- 
sen to give the lowest expected signal branching fraction 
upper limit for the null hypothesis (a true branching frac- 
tion of zero) using the MC samples. The Fisher discrim- 
inant coefflcients were determined before applying the 
|cos6'hGi|, mass, and Am requirements. We then 
tested a total of 2700 configurations of | cos6'hoi|, -7^, 
mass, and Am criteria. Table |T] summarizes the result- 
ing best values for the maximum | cos0hci|, minimum J^, 
to(Z3°) signal window, and Am interval. 

After the selection criteria in TableUwere determined, 
the data yields in the sideband region were compared 
to the expectations from Monte Carlo samples. The 
D" n+a- and D° 



e^jjL^ data yields were consis- 
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FIG. 3: Distributions of |cos(0hei)| for the three signal channels: _D° — > e^e~ (left), — >■ fi'^ fj,~ (center), and _D° — ^ e^fj,^ 
(right). The top distributions show Monte Carlo distributions for the combinatorial BB (dashed, blue) and continuum (dotted, 
red) backgrounds. The bottom distributions show the signal Monte Carlo with arbitrary normalization. 



tent with the expectations from the Monte Carlo samples. 
However, the D° e+e^ sideband yield showed a sub- 
stantial excess of events; 90 events were observed when 
5.5 ± 1.6 were expected. 

The excess of data sideband events over the expected 
background from Monte Carlo was investigated and 
found to have several distinct features: low track multi- 
plicity, continuum-like event shape characteristics, tracks 
consistent with electrons produced in photon conversions, 
low D'^ daughter track momenta, and undetected energy 
along the beam axis. We found that such events result 
from hard initial state radiation events or two-photon 
interaction processes that are not simulated in the con- 
tinuum MC samples used in the analysis. The follow- 
ing selection criteria were added in order to remove such 
background contributions: 

• Events must have at least 5 tracks for the D'-* 
e^er channel and at least 4 tracks for the 

and — >■ e^/i^ channels. 

• Events can have at most 3 electron candidates. 

• The longitudinal boost of the event, reconstructed 
from all tracks and neutral clusters, along the high- 
energy beam direction jizjE in the CM frame must 
be greater than -0.5 for all three — )■ t^l~ chan- 
nels. 

• For — > fi^iJ,^ and e^/i^ candidates, the 
pion track from the D*+ decay and the leptons 
must be inconsistent with originating from a pho- 
ton conversion. 

The signal efficiencies for the D'^ — ^ e+e^ , — > /i^/i^ , 
and D'^ — ^ e^fi^ channels for these additional criteria are 
91.4%, 99.3%, and 96.8%, respectively. The L>° ^ e+e" 
sideband yield in the data with these criteria applied is 
reduced to 8 events where 4.5 ± 1.3 are expected, based 
on the Monte Carlo samples. 



A. Peaking D° 7r~'~7r background estimation 

The amount of 13° tt+tt^ peaking background 
within the m{D'^) signal window is estimated from data 
and calculated separately for each — >■ £~^£~ channel 
using 

iVf « = (^Y. • (Ph) (P7.)) ■ 'rniDO) ■ G (1) 

where the sum i is over the six data-taking periods, N^J'^ 
is the number of D'^ tt+tt^ events that pass all of 
the -> t^i~ selection criteria except for the lepton 
identification and •m{D^) signal window requirements, 
{p'fi){pji) is the product of the average probability that 
the 7r+ and the tt~ pass the lepton identification crite- 
ria, em(_DO) is the efficiency for D'^ — >■ tt+tt" background 
to satisfy the m{D°) signal window requirement, and G 
takes into account a positive correlation in the probabil- 
ity that the 7r+ and the 7r~ pass the muon identification 
criteria. The value of (py i) {{pj i)) is measured using 
the ratio of the — tt+tt" yield requiring that the 
TT"*" (tt") satisfy the lepton identification requirements to 
the D'^ TT+TT^ yield with no lepton identification re- 
quirements applied. The {p~f j) and (pjj) are measured 
separately for each of the six major data-taking periods 
due to the changing IFR performance with time. The 
values of (Pj j) and {pj^) vary between 0.5% and 1.5%. 
The probability that the 7r+ and 7r~ both pass the muon 
identification criteria is enhanced when the two tracks 
curve toward each other, instead of away from each other, 
in the plane perpendicular to the beam axis. We use 
G = 1.19 ± 0.05 for the D° ^I+fl- channel and G = 1 
for the D° e+e~ and D° e^^i^ channels. The G 
factor is measured using a high-statistics A'^7r+ 
sample where the K is required to have a signature in the 
IFR that matches that of a tt which passes the /i identi- 
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fication criteria. This is in good agreement with the MC 
estimate of the G factor value, 1.20 ± 0.10. 



B. Combinatorial background estimation 

The combinatorial background is estimated by using 
the number of observed events in a sideband region and 
the expected ratio of events i?cb in the signal and side- 
band regions, determined from MC simulation. The side- 
band is above the signal region in the mass ([1.90, 
2.05] GeV) in a wide Am window ([0.141, 0.149] GeV). 
We fit the mass and Am projections of the combina- 
torial background MC using 2"'^-order polynomials. A 
two-dimensional probability density function (PDF) is 
formed by multiplying the one-dimensional PDFs, as- 
suming the variables are uncorrelated. The combina- 
torial background signal-to-sideband ratio Rch is then 
computed from the ratio of the integrals of the two- 
dimensional PDF. 



fit, Etttt and e^^ are the efficiencies for the corresponding 
decay modes, B^^ = (1.400 ± 0.026) x lO'^ is the 
7r+7r~ branching fraction [l^], and Su is defined by 



Su 



(3) 



The expected observed number of events in the signal 
region is given by 



iVobs = Bu/Su + Nbg- 



(4) 



III. 



RESULTS 



The uncertainties on Su and Nbg are incorporated into a 
likelihood function by convolving a Poisson PDF in iVobs 
with Gaussian PDFs in Su and Nbg- We determine 90% 
confidence level intervals using the likelihood ratio order- 
ing principle of Feldman and Cousins [T3] to construct 
the confidence belts. The estimated branching fractions 
and one standard deviation uncertainties are determined 
from the values of Bu that maximize the likelihood and 
give a change of 0.5 in the log likelihood relative to the 
maximum, respectively. 



The distribution of Am vs D'^ mass as well as pro- 
jections of Am and the D° mass for the data events for 
the three signal channels are shown in Fig. |4l Peaks 
from D'^ K~TT^ and D'^ tt+tt^ are visible at 
1.77 GeV and 1.85 GeV in the _D" mass distribution for 
— )■ fi^fJ,~ candidates. We observe 1, 8, 2 events in 
the 1?° e+e" , D" fJ.^fJ.~ , and D° e='=^'F signal 
regions, respectively. 

A. D" £'^£~ Branching fractions 

The yield of I?" — tt+tt" decays in the tttt control 
sample, selected with the same and | cos^hcil criteria 
for each Z?" — > signal mode (see Table H]), is used 

to normalize the — >■ i^£~ signal branching fraction. 
For each £>° — > i^£^ signal channel, the £>° — >• tt+tt" 
yield N^^ is determined by fitting the D'^ mass spectrum 
of the D'^ — )■ TT+TT" control sample in the range [1.7, 
2.0] GeV. The fit has three components: D'^ — tt+tt" 
, A'^TT^ , and combinatorial background. The 

PDF for the — )■ tt+tt" component is the sum of a 
Crystal Ball function and two Gaussians. The Crystal 
Ball function is a Gaussian modified to have an extended, 
power-law tail on the low side [l^l. The PDF for the 
D° — >■ Ar"TT+ component is the sum of a Crystal Ball 
function and an exponential function. The combinatorial 
background PDF is an exponential function. 

The — > £^£^ branching fraction is given by 

^« = (^) (^) - Su ■ Nu (2) 

where Nu is the number of D° £+^" signal candidates, 
N^l. is the number of tt+tt" candidates from the 



B. Systematic uncertainties 

Table HIl summarizes the systematic uncertainties. Sev- 
eral of the uncertainties in t-KixIf-ii cancel, including 
tracking efficiency for the daughters, slow pion ef- 
ficiency, and the efficiencies of the T and momentum 
requirements. The uncertainty on e^^^/ eu due to particle 
identification is 4%. Bremsstrahlung creates a low-side 
tail in the D'^ mass distributions for the 13" — > e+e" 
and — > e+/i+ decay modes. The uncertainty eu due 
to the modeling of this tail is 3% for e+e" and 

2% for — )■ e+/Lt+ . The Crystal Ball shape parameters 
that describe the low-side tail of the mass distribution 
were varied, leading to an uncertainty of 1.1% to 1.3% 
on A^* . We use the world average for the D° — >■ tt+tt" 
branching fraction [l^, which has an uncertainty of 1.9%. 
We combine the above relative uncertainties in quadra- 
ture resulting in 4.6% to 5.4% systematic uncertainties 
on Su . 

The I?" mass range for the fit used to determine the 
combinatorial background PDF was varied from [1.70, 
2.05] GeV to [1.80, 2.05] GeV. The difference in the re- 
sulting signal-to-sideband ratio Rch is taken as a system- 
atic uncertainty. The pion misidentification probabilities 
for e and /i measured in data are in good agreement with 
the MC simulation. We use the larger of either the dif- 
ference between the data and the MC or the statistical 
uncertainty on the MC misidentification probabilities as 
a systematic uncertainty. For the D° — ^ decay 
mode, we take the uncertainty on the MC estimate for 
the G factor of 8% as a systematic uncertainty on the G 
estimate from the — >■ A'"tt+ data control sample. 
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FIG. 4: Data distributions of Am vs the reconstructed _D° mass (top row) and projections of the mass (middle row) and 
Am (bottom row). The columns contain the distributions for the 7?" e^e~ (left), D" — >■ (center), and Z?" e^fj,^ 

(right) decay modes. The shaded mass (Am) distributions represent the subset of events that fall in the Am (D'^ mass) 
signal window. In the top row, the dotted (black) box indicates the signal region and the dashed (red) box indicates the 
sideband region. In the middle and bottom rows, the vertical dotted black lines indicate the boundaries of the signal region. 



C. Branching Fraction Results 

Table|TTT]prescnts the results, where Nsb is the number 
of events in the upper sideband. Neb is the expected num- 
ber of combinatorial background events in the signal virin- 
dow, is the number of events from the D° tt+tt^ 

peaking background, and Nbg (data) is the expected 
number of total background events in the data. 

For the D'^ — ?> e+e^ and D'^ — > e*/i^ channels, the 
event yield in the signal region is consistent with back- 
ground only. We observe 1 and 2 events with expected 
backgrounds of 1.0 ± 0.5 and 1.4 ± 0.3 events for the 
e+e~ and — e^ji^ channels, respectively. 
The 90% confidence interval upper limits for the branch- 



ing fractions are < 1.7 x 10 ^ for e+e and 

< 3.3 X 10"^ for D° ^ e±/^=F . 

For the — > ii^n^ channel, we observe 8 events 
in the signal region, where we expect 3.9 ± 0.6 back- 
groimd events. There is a cluster of of — 
candidate events in Fig. 2] just above and below the 
lower D'^ mass edge of the signal region, where the 
— > 7r+7r~ background is expected. We expect 7.5±0.8 
£)° TT+n- events in the entire [1.7, 2.05] GeV D° 
mass range, with 93% of these events falling within the 
narrower [1.830,1.875] GeV range. The combinatorial 
background in the [1.830,1.875] GeV D° mass interval 
is expected to be 1.8 ± 0.6 events, giving a total ex- 
pected background of 8.8 ± 1.1 events. In this interval, 



11 



TABLE II: Systematic uncertainties. The uncertainty on Su results from the uncertainties on e,r7r/e««, A^ttJt, and Btttt added in 
quadrature. The systematic uncertainty on the overaU background Nbg is obtained from the uncertainties on A'^^ and A'^^b 
added in quadrature. 





D° ^ e+e- 






j ^gg^ particle ID 


4% 


4% 


4% 


'-TTTT / ^iLk. 1 -L* J- dlioO LI Oilli 


O /O 




2% 


' TTTT 


1.2% 


1.3% 


1.1% 




1.9% 


1.9% 


1.9% 




5.4% 


4.6% 


5.0% 


^ TTTT 


11% (0.004 events) 


16% (0.43 events) 


5% (0.02 events) 


A^cb, i?cb 


36% (0.35 events) 


20% (0.25 events) 


19% (0.20 events) 


Nbg 


0.35 events 


0.50 events 


0.20 events 



TABLE III: Results for the observed event yields (A'obs), estimated background (Nbg), and signal branching fractions (Bee)- 
The first uncertainty is statistical and the second systematic. Nsb is the observed number of events in the sideband, 7?cb is 
the signal-to-sideband ratio for combinatorial background, A^cb and are the estimated combinatorial and D" — > ■k'^-k~ 

backgrounds in the signal region, A'^^^ is the fitted yield in the D° tv^tt~ control sample, e,r7r and eee are the tttt control 
sample and signal selection efficiencies, determined from Monte Carlo samples, which have negligible statistical uncertainties. 
The systematic uncertainty on e,r7r/e<f is included in the systematic uncertainty on Su, which is defined in Eqn. ([3}. 





i^o ^ e+e- 


^ /i+^i" 




Nsb 


8 


27 


24 


Rch 


0.121 ±0.023 ±0.044 


0.046 ± 0.005 ± 0.009 


0.042 ± 0.006 ± 0.008 


A^cb 


0.97 ±0.39 ±0.35 


1.24 ±0.27 ±0.25 


1.00 ±0.25 ±0.20 


^BG 


0.037 ±0.012 ±0.004 


2.64 ±0.22 ±0.43 


0.42 ± 0.08 ± 0.02 


Nbg 


1.01 ±0.39 ±0.35 


3.88 ±0.35 ±0.50 


1.42 ±0.26 ±0.20 


' TTTT 


39930 ± 210 ± 490 


51800 ± 240 ± 660 


39840 ± 210 ± 430 




14.4% 


18.7% 


14.6% 


eu 


9.48% 


6.29% 


6.97% 


Su (xlO-'') 


53.4 ±0.2 ±2.9 


80.6 ±0.4 ±3.7 


73.9 ±0.4 ±3.7 


A^obs 


1 


8 


2 


Bu (xlO-^) 


1 +"'^ 

'^■^ -0.4 


o o +2.6 
O.O _2 


5 +^-^ 

U-O -0.9 


Bu (xlO"'') 90% C.I. 


< 1.7 


[0.6,8.1] 


< 3.3 



we observe 15 events. The probability of observing 15 or 
more events when 8.8 ± 1.1 events are expected is 4.6%, 
which corresponds to a 1.7 standard deviation upward 
fluctuation from the mean for a Gaussian distribution 
(i.e. (Erf (1.7/^2) + l)/2 = 1 - 0.046). The probabil- 
ity of observing 8 events when 3.9 ± 0.6 events are ex- 
pected is 5.4%. We conclude that the excess over the 
expected background is not statistically significant. The 
Feldman-Cousins method results in a two-sided 90% con- 
fidence interval for the — > fJ.'^ fJ.~ branching fraction 
of [0.6,8.1] X IQ-'^. 

In summary, we have searched for the leptonic charm 
decays £>" e+e" , D° -)■ /Lt+^~ , and D° 6=*=^^ 
using 468 fb~^ of integrated luminosity recorded by the 



BaBAR experiment. We find no statistically significant 
excess over the expected background. These results su- 
persede our previous results [15| | and are consistent with 
the results of the Belle experiment [l^, which has set 
90% confidence level upper limits of < 0.79 x 10"'', 
< 1.4 X 10"^ and < 2.6 x 10"^ for the D° e+e" 
, — IJ-^fJ-^ , and D'^ e^fi^ branching fractions, re- 
spectively. The LHCb experiment has recently presented 
prehminary search results [l7| for — > IJ,~^ , where 
they find no evidence for this decay and set an upper 
hniit on the branching fraction of < 1.3 x 10~^ at 95% 
C.L. 
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